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q^ \ We show that a Majorana heavy neutrino with a mass 0(100 TeV) is 

Q\ | a good candidate particle for cold dark matter. It can be responsible for 

^ the majority of the cosmological diffuse photon background owing to lifetime 

Y" 1 ' of the order of (9(10 25 s), dominantly fixed by the radiative two-body decay. 

The lifetime is suppressed by two mechanisms: the leptonic GIM cancellation 

and the see-saw weak coupling suppression. As a fermion cold dark matter 

particle, a heavy neutrino favours the average mass density of the Universe 

J^j ■ constrained by the Einstein-Cartan cosmology. 



1 Introduction and motivation 

During the past decades we have been witnessed of great progress in cos- 
mology which places severe constraints on the physics of fundamental con- 
stituents of matter. There is a common belief that it is not possible to 
resolve the problems of structure formation of the Universe, basic cosmolog- 
ical parameters, nucleosynthesis or a diffuse photon background without the 
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introduction of cold and hot dark matter and the violation of the baryon 
number. 

On the other hand, measurements of the LSND and the SuperKamiokande 
have definitely confirmed the existence of massive neutrinos and their flavour 
mixing, urging our necessity to build a more predictable theory in particle 
physics than the Standard Model. 

A diffuse photon background as a source of ionization at galactic or ex- 
tragalactic scales represents a challenge for particle physics to search for a 
process that can explain the phenomenon |Q. It has been known for a long 
time that the measurement of the flux of decay-produced photons can con- 
strain the lifetime of the decayed particle 0, || |J. From a detailed study of 
ionization fluxes |I||, Sciama has estimated that the mass and the lifetime of 
a decaying light neutrino are m„ T — 28eV, r ~ 10 24 s. However, it is very dif- 
ficult to reconcile such a high mass with the present cosmological fits of the 
structure formation which require an order of magnitude smaller neutrino 
masses (as hot dark matter) and cold dark matter (that should dominate 
the mass density of the Universe). The observed neutrino oscillations at the 
SuperKamiokande and the limits from other detectors prohibit large masses 
if we exclude mass degeneracy. The standard weak interaction calculation of 
the neutrino lifetime|| with a mass m v ~ "5§eV gives r ~ 10 36 s, thus twelve 
orders of magnitude larger than it is required by the ionization flux. One can 
improve the result adding a contribution of a new scalar or vector charged 
particle exchanged in the quantum loop, but new particles could completely 
spoil the structure of electroweak interactions that are verified to very high 
precision at LEP and SLD. To measure a photon flux of the decaying neu- 
trino near the Sun, an extreme ultraviolet detector on the satellite has been 
proposed 0. In 1997 the satellite was launched and now there are data that 
are incompatible with Sciama's decaying neutrino ||. 

It seems that light neutrinos are not very promising candidates to solve 
the cosmological problem of diffuse photon background. 

In this paper we investigate abundances and lifetimes of heavy neutrinos, 
within an electroweak theory proposed a few years ago, in order to find a 
cold dark matter(CDM) particle whose radiative decay could be responsible 
for a diffuse photon background. 



2 Cosmic abundances of heavy neutrinos 

To investigate the cosmological significance of any particle, one should know 
its interactions and cross sections to solve the Boltzmann equations in curved 
spacetime. Let us recall some results on the freeze-out and abundances of 
heavy neutrinos. 

Heavy neutrinos with masses 0(lGeV) can play a role of the CDM par- 
ticle owing to high abundances calculated from the Z-boson mediated anni- 
hilation cross section [[J with the following scaling oz oc m 2 N . If neutrinos 
have masses much larger than weak bosons, the authors of Ref. showed 
that the total annihilation cross section was not dominated by the fermion 
pair production a(NN — > //) oc m^ 2 , but by the W-boson pair production 
a(NN — > W~W + ) oc m 2 N . They concluded || that "there is no cosmological 
upper bound on the masses of very heavy stable neutrinos" . They assumed 
the standard weak coupling of (Dirac) neutrinos. 

However, Griest and Kamionkowski successfully put an upper bound 
on the mass of the CDM particle estimating the upper bound of the total an- 
nihilation cross section without reference to any particular interaction model 
or theory. 

Nevertheless, it should be remembered that heavy neutrinos have a sub- 
stantially different scaling of the cross sections in the fermion or boson pair 
production. If one excludes the W~W + production, one can reach a cosmo- 
logically acceptable abundance of a heavy neutrino as a dark matter particle 
with a mass m N ~ 900GeV [|]. 

It is now necessary to define our framework for a study of heavy neutrinos. 
We work within an electroweak theory proposed some years ago H|. We 
called it (BY) theory in that paper. It differs from the Standard Model(SM) 
(called (AX) theory) in two essential ingredientes: (1) Nambu-Goldstone 
scalars carry nonvanishing lepton numbers and only Majorana fields acquire 
masses at tree level, (2) there is no Higgs scalar and the principle of the 
noncontractible space (the existence of finite UV scale) is introduced into the 
local gauge field theory as a symmetry-breaking mechanism, fixing the tree- 
level weak boson and Majorana fermion masses. The theory contains three 
light and three heavy Majorana neutrinos. Further insights and discussions 
could be found in Ref. 0. Here we concentrate on the leptonic sector. 

Let us recall the structure of the four-component Dirac spinor: 



*!> = (£> (i) 



n a transforms under a matrix S of SL(2,C), 



£ a transforms under a matrix (S )* of SL(2,C). 

To retain all degrees of freedom, one can choose the following set of 
fermion fields: 



e L , e R , * L , (V c )l. (2) 

These fields interact with 577(2) x [7(1) gauge fields and Nambu-Goldstone 
scalars as follows: 

£ — £-lep + £-g.bos + £-scal + ^Yuk + £g.fix + '-'FPi (3) 

C lep = RiYidp + ig'B^R + Zz 7 "(0„ + \g'B^ - ig^A^L + ^^d^ R , 

r , -_IP F W _ I R RM 1 ' 

'-'g.bos 4- 1 ^i/- 1 ± ±J iw ± -' i 

Cscax = (d^ + i^B^ + zfr i ^$t)(a^$ - z^S^ - ifrM^$), 
definitions : e = charged lepton, ip = neutral Dirac leptc 

F% = dpAl - d v A\ + ge^A^Al B, u = d,B u - d u B„ 



ton, 



b°(x 



= (v + ix°(x))/v2, v = symmetry breaking parameter, 
~,X° = Nambu — Goldstone scalars. 
The physical spectrum of neutral leptons looks like |7|, § 



m L m D \ I a \ I a 

m D ) \b ~ A U 



a z + b' = l, 



^n = a N f + b N F, ■$„, =a v ,f + b v >F, (4) 



/ = ^l + (^) c ), F = -^(tt* + (* fl ) c ), 



Atv = - (m L + sjm\ + Ami), V = ^( m L - \jm\ + 4m?,), (5) 
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&„/ ~ -1, 



Oat ^ 1, b N ~ 

m L m L 



*„ = 75*^' =► A„ = -A,/. (6) 

Heavy Majorana particles acquire masses predominantly from the quan- 
tum loops with a Nambu-Goldstone scalar in the strong coupling regime of 
Dyson-Schwinger equations [^]. However, the effective strong coupling for 
fermions with masses of a few TeV is saturated in bootstrap equations and 
its value (like the QCD strong coupling in the infrared domain) is not so high. 



The cross section for heavy neutrinos annihilating into a heavy neutrino pair, 
calculated in the 't Hooft-Feynman gauge is 



a^iNiNi -► NjNj) = . 1 a {d e )i(d e )j (7) 

sm t>w 

1 s — 2m 2 N , 



Notice that we cannot decouple Nambu-Goldstone scalars when they are 
in the strong coupling regime with fermions. Only a complete solution of 
the whole set of Dyson-Schwinger equations can give us gauge-invariant ob- 
servables. Thus, the preceding cross section formula should be read off as 
only for effective quantities. Because of the see-saw suppression, the cross 
sections mediated through electroweak gauge bosons are negligible and, in 
addition, there is no W~W + pair production via the x° Nambu-Goldstone 
boson because Vertex(x°W~W + ) = and there is no Higgs scalar in the 
theory j7|. The scaling of the cross section in Eq.(7) is the same as that of 
the cross section az(NN — > //) (Eq.(2) of Ref. ||): a x o oc mjj 2 . 

Knowing the predominant contribution to the total annihilation cross 
section, we can estimate the abundances of heavy neutrinos. The freeze-out 
temperature of the CDM particle depends on the cross section only logarith- 
mically [0], so it is not neccessary to solve the freeze-out condition M for a x o. 
We can compare cross sections and make the following estimates: 



s > m 2 N] ,m x o,r x o, 



vr ,_,,_.„ 1 



(J v o 



f-\> 



(aeMOe)^-, 



sm t>w s 



(3 N = (1 - Am 2 N /s) l '\ s~4m 2 N + 6m N T f , T f 
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x f = ^, x f ~ 17 (mjv ~ 2GeV), x f ~ 25 (m N > ITeV), 



(a e ) e ~ 0.4, (a e ) M ~ 1, (a e ) T ~ 6, 



fi/v oc 



(n + l)x^ +1 

J - — , {a\v\} = a x 

ctq 



<?z(mN — 0.9TeV) ~ cr x o(mjvJ, ^jv — 1 



=► mjvi = O(lOOTey), i = fi,r. (8) 

Although we are limited with the knowledge of the effective coupling 
of fermions to the Nambu-Goldstone scalars, one can conclude that heavy 
Majorana neutrinos with a mass of order O(100TeV) are good candidates 
for the CDM particle. 

If one assumes that the radiative decay of the CDM particle dominates, 
we can estimate its lifetime from the measurements of the differential energy 
flux of the diffuse photon background || : 



n N ~ 1 and Fig. 7 of Ref.[9] 



^r N = C(10 25 s) for m N = O(100TeV). (9) 

The present estimates of the mass and lifetime of the CDM particle are 
further challenges for a theory, so the next section is devoted to a study of 
its lifetime. 



3 Lifetime of heavy neutrinos 

Searching for the tree-level decay processes of neutrinos, one can recall that 
the tree-level decays of light neutrinos are kinematically forbidden. Heavy 
neutrinos couple to electroweak gauge bosons with see-saw suppression fac- 
tors, but in our theory these factors vanish at tree level [0] because 
rriDi (tree level) = 0, i=flavour qauntum number. Just like for light neutrinos^ 
we have to calculate the two-body radiative flavour changing decay of heavy 
neutrinos at the quantum loop order. 
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Let us be concerned with the Ni — > v^ decay via the charged weak 
current loop. When dealing with processes containing light and heavy neu- 
trinos, one has to symmetrize fermion states with respect to the following 
interchange of fields: ipi ^ 4>i => fi *~* Fi> i — flavour. Nature should not 
recognize what we define as a particle or an antiparticle. Acknowledging the 
see-saw mixing coefficientes from Eq.(4) and after the particle-antiparticle 
symmetrization, one can conclude 

T(Ni -> Ujj) oc {a N .a v >. + K>b Nj ) 2 , V i,j, 

a Ni = -b v i, and a v >, = b Nj ^T(Ni -> Ujj) = 0. (10) 

We shall now study the decay processes Ni — > Nj-f induced by the loop 
exchange of the W weak boson and charged leptons. The effective transition 
operator will be estimated preturbatively, so one has to choose the renormal- 
ization conditions for the flavour mixing of leptons. 

The usual wisdom for fermion mixing is to apply the on-shell fermion 
mixing scheme of Ref . |H| . It has been shown that in flavour- changing lepton 



radiative processes, this scheme causes the appearance of the dimension-four 



operators [11] besides the standard dimension five-operators M. However, 



the authors of Ref.|12[ have recently shown that the on-shell renormaliza- 
tion scheme is not consistent because it violates Ward-Takahashi identities 
and leads to gauge-dependent physical amplitudes. Instead, one can in- 
troduce a natural and consistent prescription where flavour- changing self- 
energies vanish at zero momentum |T3]. From the Ward-Takahashi identity 



and the conservation of the electromagnetic current one can conclude that 



new renormalization conditions of Ref. |L2J do not induce the appearance of 
four- dimensional operators in flavour changing lepton radiative decays. 

The effective see-saw suppressed transition operator for iVj — > Njj, (i, j = 
flavours) is evaluated in the t'Hooft-Feynman gauge (loops with exchanged 
WWl and Wll, l=charged lepton) because only a strongly coupled heavy 
neutrino-Nambu-Goldstone scalar (decoupled in the unitary gauge) system 
can produce the mass splitting of heavy neutrinos (CP invariance assumed, 

I): 

A, = [H(0) + H 5 (0) l5 ]a, p (p 2 r, (11) 



opposite CP eigenvalues of neutrinos : 



eg* / m Vi m Vj 
167T 2 V m Ni m Nj ^p 



H = iL J^E^ E UtMHL + H R ),H 5 = 0, 



same CP eigenvalues of neutrinos : 



h 5 = ^ 2 J^^j:u; j u u (h r -h l ),h = o, 

\m A y m Ni m Nj ~ J 

H L = m Nj (C + C n + C 12 + C 23 ), 

Hr = m Ni (—C + C12 + C12 — 2Cn — -Cn — C 2 i + C23), 
Co = C (pl,pl,pl;M w ,mi,mi), 

C = C (pl,i%,pl;rrn,M w ,Mw), 

p\ = m 2 N ., p\ = 0, p\ = m 2 N .,m[ = charged lepton mass, 
for further definitions see Appendix A. 

The amplitude is ultraviolet (UV) finite, but it contains infrared(IR) sin- 
gularity in the limes mi — > 0. It can be visualized that the IR singularity 
comes from the following two Green's functions: 



3tB (0;mi,mi) = — lnm^ + 
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$tC = -5 K(-ln 2 wi - -ln 2 w 3 ) + ..., 

p( -pi 2 2 



w(p 2 ,M w , mi ) = — (p 2 -M 2 w + m 2 - yV - K? + mf) 2 - 4p 2 m 2 ), 
zp 

Wl = w(p 2 ,M w ,mi), w 3 = w(pl,M w ,m l ). 

The IR singularity can be removed in a natural way through the leptonic 
GIM mechanism: 



H = Hm £[W(mf, ...) - H(m 2 = M 2 , ...)], 

similarly for H 5 . 



The IR singular and constant terms of the whole amplitude are now 
subtracted away. 

Instead of the Higgs mechanism, the finite UV scale is introduced into 
the theory, so one has to study the finite-scale effects. They enter into the 
calculations through the evaluations of scalar Green's functions |H| . One can 
naively expect large corrections, but the explicit evaluation (see Appendix B) 
tells us that they are three orders of magnitude smaller, thus much smaller 
than uncertainties in masses. We shall neglect this effect in the numerical 
estimates of lifetimes. The reason for small corrections is that for these decay 
processes we need the knowledge of the Green's functions in the timelike 
region for very high masses (m^ ^> A). The contribution from the timelike 
region dominates over that of the spacelike region. On the contrary, in QCD 
one studies the spacelike region with a cutoff A, thus a large deviation will 
be encountered [TBI for the scale fi > A. 



A straightforward calculation gives the partial decay width of the Majo- 
rana heavy neutrino (assumed the same CP eigenvalues of neutrinos ||): 



1777 777 1 

T ( Ni -> N j7 ) = ^ % 3 ^ | H 5 | 2 . (12) 

8irm 3 N . 

Because of the fact that m 2 N . 3> M^, a partial decay width is insensitive 
to the weak boson mass to the leading order: 



T(N t -> Nfy) oc m Ni sin 2 (2%)(^^)(^) 4 . (13) 

m Ni m Nj m N] 

Owing to this scaling, the gauge dependence of the effective transition 
operator is of subleading order and there is no need to study the gauge can- 
cellations in detail. The heavy-quark symmetry of Isgur and Wise or the 
Appelquist-Carazzone decoupling theorem are also examples in field theory 
where a good approximation in a certain sector of the complete theory de- 
scribes the actual physical situation correctly. 

Next section is devoted to final numerical evaluations and physical con- 
clusions. 
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4 Results and conclusions 

From the requirement for the abundance of the CDM particle and the require- 
ment of the cosmic diffuse photon background we extracted the information 
about possible masses and lifetimes of heavy neutrinos. The general scal- 
ing of the lifetime on light- and heavy-neutrino massees and mixing angles 
suggests that there is no dependence on the mass of the decaying particle, 
because of the cancellation(see Eq.(13)). 

Now we present our numerical results with an appropriate choice of neu- 
trino masses and mixing angles. From the scaling relation Eq.(13) one can 
easily estimate a partial decay width for some different set of parameters. 



m Ne = lOTeV, m N = 50TeV, m Nr = 200TeV, 



m Ve = 0.03eV, m v = 0.3eV, m VT = 2eV, 



012 = 0.21, 23 = 0.02, 31 = 0.002, 



Er- 1 ]- 1 , 



t(N t -> JV 7) = 2.3 • 10 25 s, t(N t -► N el ) = 5.4 • 10 24 s, 



t(N^ -> N el ) = 2.0 ■ 10 26 s. (14) 

The lack of our knowledge of the mixing angles and neutrino masses does 
not prevent us from drawing main conclusions of this paper concerning heavy 
neutrinos: 

(1) Acknowledging the see-saw relation for light-neutrino masses and the re- 
lation for the mixing angles 9\y = 2(6>i 2 + 023 + 03i) 0, the (BY) theory 
predicts that the heavy neutrino N^ or iV T with the mass O(100TeV) and 
the lifetime r = O(10 25 )s is a CDM particle which can solve the problem 
of the cosmic diffuse photon background; this means that a further detailed 
astronomical study of the diffuse photon background could be useful for the 
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study of the dynamics of the CDM particle, which could be supplemented by 
terrestrial measurements (for example, the DAMA experiment at the Gran 
Sasso National Laboratory); it is not excluded that processes with heavy 
neutrinos could explain the recently observed(EGRET) diffuse gamma halo 
around our Galaxy. 

(2) Essential ingredients of the (BY) theory, which gives successful phe- 
nomenological answers, are the peculiar chiral structure of the theory, the 
absence of the Higgs scalar and the presence of the finite UV scale to fix 
the scale of heavy neutrinos and weak gauge bosons, as well as the presence 
of the UV finite self-consistent bootstrap system of equations leading to the 
finite number of fermion families. 

(3) Because of the rather strong tree-level coupling of N% to x°'- j^fg (j^) 2 = 
0.29, m = 485GeV, it seems that even the lightest one of heavy neutrinos 
acquires the mass mostly from the quantum loop (m^ e 3> mo). 

(4) The forthcoming higher-luminosity measurements at Tevatron could con- 
firm the nonresonant enhancement of the QCD amplitudes beginning in the 
vicinity of A (lim^oo af (//) ^ 0) [[L3|]. Also note the unsettled problem of the 
anomalous b- quark weak coupling (LEP and SLD) as a possible consequence 
of the finite scale effect in the loop corrections with the t-quark exchange 

(5) The Einstein-Cartan cosmology requires that Q m ~ 2, thus the CDM 
particle must be the spinning particle (fulfilled for the (BY) theory) |L4|] ; 
this setting of the average mass density of the Universe is a consequence 
of the following bootstrap at the level of the Einstein-Cartan equations: the 
number density of the particle that dominates the energy momentum density, 
at the same time dominates the spin density of the Universe, but the energy 
momentum and spin of matter are coupled to the curvature and torsion of 
spacetime with the same coupling constant. 



(6) The current controversies [15] concerning the average mass density of the 
Universe could be resolved only with the most general cosmological model 
that contains cold and hot dark matter, cosmological constant, baryons and 
CMBR and the metric not only with expansion, but also with a small amount 
of vorticity, acceleration and shear f[4jl ; namely, one can solve the primor- 



dial mass density fluctuation within the Einstein-Cartan cosmology only be- 
yond the standard Friedmann-Lemaitre- Robertson- Walker metric ||16| : it was 



shown a long time ago that the primordial vorticity could generate a cosmic 



magnetic field |T7 
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Appendix A 

Here we display the scalar and tensor Green's functions used in the evaluation 
of the effective transition flavour-changing operator: 



A(m) = —l'd A q- ' 



m 2 J q 2 — m 2 + is ' 

B (p 2 ;m 1 ,m 2 ) = — » / A71 —, ^77— — r» ~, V' 

m z J {q A — m{ + t£){{q + p) — m 2 + is) 

B (v 2 -m m ) = — — I d A a — 

M ' ~ nr 2 J (q 2 — m 2 + ie)((q + p) 2 — m 2 + ieY 

g^B 2 2 +P^PvB 2 l = 2 / ^ Q~T2 2~, \(( 1 7 \2 2~, 7' 

itx a J [q z — m{ + i£){{q + p) — vn 2 + is) 

Pl+P2+P3 = 0, 

d A q 



C (p 2 1 ,pl,pl;m 1 ,m 2 ,m 3 ) = — \ 

ITT J 



(q 2 — m\ + ie)((q + Pi) 2 — m 2 + ie)((q — P3) 2 — rn 2 + is) ' 
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PiCn + p 2 C 12 = —3 / Ati 2" 

271^ j (g^ — mf 



+ «e)((g + pi) 2 -ml + ie){{q - p 3 ) 2 + m\ + te) 



gTCn + p^C 21 + p^C^ + (pfr£ + p^)C 23 



ITT" 



d 4 q 



qiiqv 



(q 2 — m 2 + ie){{q + Pi) 2 — m 2 + ie)((q — p 3 ) 2 + m 3 + is) ' 



B x {p ;m 1 ,m 2 ) = 7 —^[A(m 1 ) - A(m 2 ) + (m 2 -m x -p )B (p ;m 1 ,m 2 )} 



2p 



1 



B 22 (p 2 ; mi ,m 2 ) = ^[A(m 2 ) + 2m\B Q {p 2 ; m u m 2 ) + (p 2 + m\- m 2 2 )B l {p 2 - m u ma)], 

B 21 (p 2 ;m 1 ,m 2 ) = — -[A(m 2 ) - m 2 1 B (p 2 ;m 1 ,m 2 ) -2(p 2 + m\ - ml)B 1 (p 2 ;m 1 ,m 2 )}, 
ip* 



C 



12 



(Pi ■ P2) 2 ~ P1P2 



\Pl^2 + (pl + Pi -P2)S 



1J- 



C U = -o(Ei -pi -P2C12) 
Pl 
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Si = -B (p 2 3 ;m 1 ,m 3 ) + -^{~p\ ~ m\ + m 2 2 )C (p 2 1 ,p 2 2 ,p 2 3 ;m u m 2 ,m 3 ) - -B (p 2 2 ; m 2 , m 3 ) , 
S 2 = -iBb(p?; mx, m 2 ) + ±(pg - m\ + ro?)C (p?,jfc pg; m b m 2 , m 3 ) + ±S (p& m 2 , m,), 



C, 



2:; 



2(PiPl - (Pi • P2)' 



■[(-Spl - 4pi ■ p 2 )A 2 - 3p?A 4 + Pl ■ p 2 A ± - p 2 2 A 3 



C, 



1 



21 



P\ -Pi 



[A2 + A4+P2C23], 



C 22 = [A 3 -p\C 23 } 

P\-P2 
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C 2 4 = A 2 - plC 2 i - pi • p 2 C 2 3, 
Ai = 5 (pi; ma, m 3 ) + m 2 C , 



A. = i Bl «; mi , m .) + i(-ri + ^- m SC I1 , 



A 3 = gSi^; mi, m 3 ) - -B 1 (pl;m 2 ,m 3 ) + -(-pj + mj - ml)C 12 , 



Appendix B 

The real parts of the two- and three-point scalar Green's functions in the 
noncontractible space are given as in Ref. fll3|: 



KE A (p 2 ; mi, ma) = ( / c^(p 2 , j/) + % 2 - m 2 ) / ^AK(p 2 , y)) 

JO J — (\/p -112) 

^(p 2 ,y) = 



(Vp 2 -m 2 )2 ' ' y + mf 

2y 



p 2 + y + m 2 . + \/{—p 2 + y + m 2 ,) 2 + 4p 2 y 



Air(p 2 ,?/) 



-p 2 + y + m 2 ,) 2 + 4p 2 y 
p 2 



The integration in the second term is performed from the branch point of 



the square root \/(—p 2 + y + m 2 .) 2 + 4p 2 y = iZ and the additional kernel is 

-p 2 +y+rri2+iZ 



derived as the difference: Ai^(p 2 , y) = K(p 2 , y) — K*(p 2 , y) — — 



Ja_ 



— p 2 +y+m2— iZ ' 

The integration over singularities is supposed to be the principal-value 
integration. 

In the case of the two-point Green's function Bq, we need the explicit 
form of the additional term for the integration in the timelike region because 
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the integration in the spacelike region is divergent in the limes A — > oo. 
However, the three-point scalar Green's functions are UV-convergent and we 
do not need to know the explicit form of the additional terms because they do 
not depend on the UV cut-off and we can use the analytical continuation of 



the standard Green's functions written in terms of the dilogarithms [ 13| |18 



^CQ(p i} m j ) = dq 2 ^(q 2 ,p i ,m j ) + dq 2 E(q 2 ,p i ,m j ), 

Jo Jtd 

3tC£(pi, rrij) = ReC^ipi, rrij) - / dq 2 $(q 2 ,Pi, rrij), 

J A 2 

$ = function derived by the angular integration after Wick's rotation, 
C^ = standard 't Hooft — Veltman scalar function, 
TD = timelike domain of integration. 

This equation is valid for arbitrary external momenta. The same formula 
is applicable to the higher n-point one loop scalar Green's functions. 

To confirm the claim that the finite scale effects in the heavy-neutrino 
radiative decay are small, we evalute the C = Co Green's function as an 
example: 



^C A = $tC°° - AT, 



Ar = r°° - r A , 



2 f A , q 



n Jo q 2 + M^ 



2; 



IV 



r+1 

I\ = I dx 



-mi x 



l (q 2 + m 2 + mf) 2 + 4m 2 q 2 x 2 k 



q 2 + ml + mf + 2kq^/l — x 2 q + ml + mf — 2kq^/l — x 2 . 

x (arctan — — : : arctan — — : : 

2(mi + k)q \ x \ 2(mi + k)q \ x 
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f +1 1 q 2 + m\ + mf 

1 (q 2 + m\ + mf) 2 + Am\q 2 x 2 Akq 



(q 2 + m\ + mf + 2kq\J\ — x 2 ) 2 + 4(mi + k) 2 q 2 x 2 
(q 2 + m 2 + mf — 2kqy/l — x 2 ) 2 + 4(mi + k) 2 q 2 x 2 



mi — ml 

k = — -, 

2m\ 



pl = m \ = (lOTeV) 2 , v \ = m\ = (0.5TeV) 2 , A = 326GeV, 
AKC°° = $tC°°{mi = m„) - $tC°°(mi = m e ) = 1.00 • lO^TeV 
A(Ar) = AT( mi = m M ) - AT( mi = m e ) = -6.87 ■ 10~ 12 Tel/~ 

=H A5RC°° |>| A(Ar) | . 
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